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We introduce and solve analytically a model for the development of disparate social classes in 
a competitive population. Individuals advance their fitness by competing against those in lower 
classes, and in parallel, individuals decline due to inactivity. We find a phase transition from a 
homogeneous, single-class society to a hierarchical, multi-class society. In the latter case, a finite 
fraction of the population still belongs to the lower class, and the rest of the population is in the 
middle class, on top of which lies a tiny upper class. While the lower class is static and poor, the 
middle class is upwardly mobile. 
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A distinguishing feature of developed societies is the 
existence of social classes. What accounts for this diver- 
sity? Do individuals advance their position as a result of 
innate talent, inherited wealth, or plain luck? Our aim 
is to construct a minimalist interacting agent model that 
accounts for the development of social diversity. 

The phenomenon of social diversity has been exten- 
sively investigated both in the social 0, Q] and the bi- 
ological sciences 0, • Social hierarchies have been ob- 
served empirically in a wide range of animal populations, 
from insects to mammals 0, and to primates and 
humans An appealing route for modelin g so cial dy- 
namics is to use interacting agent systems M llOj . Quan- 
titative an d large scale simulations |l2| of physics- 
inspired interacting particle systems have been used to 
model observed macroscopic collective phenomena using 
microscopic agent interactions [ljj llj, ll^i • 

Motivated by empirical observations of bumblebee 
communities pjj, Bonabeau et al. rece ntly introduced 
an agent-based model of social diversity [17J where each 
individual is endowed with a fitness-like variable that 
evolves by two opposing processes. The first is competi- 
tion: when two agents interact, one individual becomes 
more fit (gains status) and the other becomes less fit, with 
the initially fitter individual being more likely to win. 
Counterbalancing the competition, the winning probabil- 
ity for the fitter agent decreases as the time from the last 
competition increases. Investigations of this model have 
found a transition from a homogeneous to a hierarchical 
society as the relative strengths of these two processes 
are varied [HEEi. 

In this letter, we introduce a simplified version of the 
Bonabeau social diversity model that accounts for the 
competing processes of advancement by competition and 
decline by inactivity via a single parameter. By solving 
the underlying rate equations, we find a phase transition 
from a homogeneous, single-class society to a hierarchi- 
cal, multi-class society. In the latter phase, the lower 
class is destitute and static while the middle class is dy- 
namic and has a continuous upward mobility. 

In our model, an agent is endowed with an integer 



fitness value k > 0. All agents start with zero fitness 
and the fitness changes by two processes: (i) advance- 
ment by competition and (ii) decline by inactivity. In 
the competition step, when two agents interact, their fit- 
ness changes according to 



(fc,j)-> (k+l,j), 



(1) 



for k > j. When two equally fit agents compete, only one 
advances. Without loss of generality, the rate of this pro- 
cess is set to one. We also consider the mean-field limit 
where any pair of agents is equally likely to interact. The 
rationale behind this "rich gets richer" dynamics is ob- 
vious: fitter individuals are better suited for, and hence 
benefit from, competition. When decline occurs, individ- 
ual fitness decreases as 



1 



(2) 



with a rate r. This process reflects the natural tendency 
for social status to decrease in the absence of interactions. 
The lower limit for the fitness is k = 0; once an individ- 
ual reaches zero fitness, there is no further decline. The 
model is characterized by a single parameter, the rate of 
decline r. 

Let fkif) be the fraction of agents with fitness k at time 
t. This distribution obeys the nonlinear master equation 



dfk 
dt 



r (/fe+l — fk) + fk-lFk-l — fkFk 



(3) 



for k > 0, where Fk = Y)j—Q fj is the cumulative 
distribution. The boundary condition is /_ i = so 
that dfo/dt = rfi — /q, and the initial condition is 



fk(0) 



The first two terms in Eq. @ account for 



the decline of the fitness of an agent, while the last two 
terms account for advancement |20|. 

To understand the behavior of this system, we focus on 
the cumulative distribution Fk , from which the individual 
densities are fk — Fk — -Ffc-i- From the master equation 
P|. the cumulative distribution satisfies 



dF k 
dt 



r{F k+ i - F k ) + F k {Fk-! - F k ), 



(4) 
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FIG. 1: The middle class. The scaled cumulative distribution 
$0) versus x for r = 1/2 at t = 250 (dotted), 1000 (dashed), 
4000 (dot-dashed). The solid line is the theoretical predic- 
tion JSJ. The inset shows the qualitative behavior for r — 
(dashed), r « 1/2 (solid), and rwl (dotted). 



for fc > 0. The boundary condition is F_i = so that 

dF /dt = r(Fi — F ) — F§, and the initial condition is has a Gaussian decay. 
F k (0) = 1. 

Homogeneous vs. Hierarchical Societies. Our so- 
cial diversity model undergoes a phase transition from 
a homogeneous to a hierarchical society. This transition 
follows from the continuum limit of the master equation 
Q for the cumulative distribution 



FIG. 2: The lower, middle, and upper classes. The fitness 
distribution /j. versus k for r — 1/2 at t = 4000 (solid line), 
showing the lower (k < 60), middle, and upper (k > 2000) 
classes. Also shown for reference is the plateau = 1/t 
(dashed line). The distribution decays as 1/k 2 in the lower 
class up to the diffusive scale fci owcr ~ (2rf) 1//2 . The distri- 
bution is constant in the middle class up to a ballistic scale 



(1 — r)t, beyond which there is an upper class that 



Eq. Q suggests that the characteristic fitness of the mid- 
dle class increases linearly with time, k ~ t. Thus, we 
posit the scaling form 



F k ~ 



(7) 



dt y ' dk 



(•5) 



For finite fitness, the cumulative distribution approaches 
a steady state in the long-time limit. Then either F = r 
or dF/dk = 0. Invoking the bound F < 1, we conclude 
that either F = r or F = 1. Therefore, L, the frac- 
tion of the population with finite fitness exhibits a phase 
transition 



with the boundary condition <f>(oo) = 1. Substi- 
tuting Eq. into (J5J, the scaling function satisfies 
xdQ/dx = ($ — r) d&/dx where x — k/t. The solution 
is either = r + x or d&/dx — 0. As a result (Fig.0 



$(x) = 



X < 1 

£ > 1 



r; 
r. 



(8) 




(6) 



When competition is weak, the entire population has a 
finite fitness, while for strong competition, only a fraction 
L < 1 of the population has a finite fitness. 

We shall see that the quantity L is the size of the lower 
class, while the complementary fraction 1 — L is the size 
of the middle class, whose fitness increases indefinitely 
Thus for r > 1, the society is homogeneous and consists 
of a single lower class. However for r < 1, there is a 
hierarchical society that contains a distinct lower class, 
and a distinct a middle class. When r = 0, the lower 
class disappears entirely. 

Middle Class Dynamics. The picture presented above 
is confirmed by analyzing the dynamics of the middle 
class. Applying dimensional analysis to the governing 



Remarkably, the scaling function for the cumulative dis- 
tribution is piecewise linear and thus non- analytic. 

The scaling function iJSJ has a number of basic implica- 
tions. First, the quantity $(0) = r is the fraction of the 
population that belongs to the lower class, confirming 
the prediction of Eq. ©. This behavior is reminiscent 
of a physical Bose condensate, where a finite fraction of 
the population occupies the zero fitness (in scaled units) 
ground state. In this sense, the entire lower class is desti- 
tute. When there is only competition (r = 0), the society 
consists of a continuously-improving middle class. In this 
case, a formal exact solution of the master equations is 
possible |2l| . 

We can alternatively write the fitness distribution in 
the scaling form fk — t _1 ^(fc/t). The corresponding scal- 
ing function is <j)(x) = d$/dx = rS(x) + 1 for x < 1 — r 
and 4>{x) — otherwise. The middle class thus has a 
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FIG. 3: The lower class. The cumulative distribution, nor- 
malized by the rate r, Fk/r is plotted versus k. Shown are 
simulation results for r = 1/4 (circles), r — 1/2 (squares), 
r = 3/4 (diamonds) at time t = 10 4 . 




FIG. 4: The upper class, 
fitness distribution: tf k 
z = (k — vt) I yDt, for r - 
t = 8000 (squares). 



Shown is the normalized tail of the 
versus z 2 , with the scaling variable 
= 1/2 at times t = 4000 (circles) and 



constant fitness distribution 



(9) 



for k < fcuppcr = (1 — r)t. The lot of the middle class is 
constantly improving, as the fitness extends over a grow- 
ing range and the average fitness increases linearly with 
time. 

Numerical integration of the master equation confirms 
these predictions (Figs.Q]and[21). We used a fourth-order 
Adams-Bashforth method |22( with accuracy to 10~ 10 in 
the distribution F k . Our numerical data was obtained by 
integrating F k for < k < 20000. 

Lower Class Dynamics. The fitness of the lower class 
is finite; in other words, the fitness distribution is in a 
steady state. This distribution can be determined by 
setting the time derivative in the rate equation to zero. 
Writing Fk = L(l — Gk), so that the deviation Gk van- 
ishes at large k, Eq. (QJ gives 



G. 



k+l 



Gk 



Gk — G 



fc-i 



L(l-G k ). 



(10) 



The fitness distribution is fundamentally different in 
the two phases. In the homogeneous society phase (r > 1 
and L = 1), the deviation Gk decays rapidly at large 
fitness. Replacing the right-hand side of Eq. (|10ll by 1 



for large k, the solution is simply G k 



.-k 



fk 



,-k 



Therefore 



(11) 



The fitness distribution decays exponentially, so that the 
lower class is confined to a small range of fitness values. 
The characteristic fitness 1 / In r diverges as the transition 
is approached. The society is homogeneous as it contains 
a single social class, the lower class, that does not evolve 
with time. 



In the hierarchical society phase, (where r < 1 and 
L = r), the fitness distribution is universal, as 
the recursion relation I jlOp becomes r-independent, 
(G k+ i - G k )/{G k - Gk-x) = 1 - G h . This shows that 
Fk/r is a universal, r-independent distribution (Fig. EJ. 
We start by treating k as a continuous variable, because 
the fitness range becomes large as r J, 1. We thus ex- 
pand the differences in Eq. (|10|l to second order. Since 
G" <C G', where prime denotes differentiation with re- 
spect to k, we find G" + GG' = 0. Integrating once and 
invoking G — > as k — > oo, gives G' + \G 2 — 0. Asymp- 
totically, G ~ 2k~ 1 , and then by using f k — F k ~ Fk-i, 
we find 



fk^2rk~ 



(12) 



The lower class has a power-law fitness distribution with 
mean fitness that diverges logarithmically in the upper 
limit. While the lower class is still static, it is not as 
destitute as in the homogeneous society phase. 

The transition between the lower and middle class oc- 
curs when 2r/k 2 l/t, i.e., where the power-law distri- 
bution H12(l matches the uniform distribution . Conse- 
quently, the lower class is confined to a diffusive boundary 
layer of thickness 



Mower 



(2rt) 



1/2 



(13) 



Beyond this diffusive scale, lies the middle class 
whose constant density extends over the range 
Slower < k < fcuppcr- In the hierarchical society phase, the 
fitness distribution consists of the stationary component 
(|12|l that defines the lower class and the evolving compo- 
nent (JZJ) that defines the middle class. The extent of the 
stationary region indefinitely grows with time. 

We thus conclude that the lower class is always static, 
being in a steady-state independent of the rate of de- 



4 



cline r. In a homogeneous society, the lower class has 
an exponentially decaying fitness distribution that lies 
within a narrow fitness range. In a hierarchical society, 
the lower class fitness distribution decays algebraically 
and its range grows diffusively with time. 

Upper Class Dynamics. The upper class is de- 
fined by the subpopulation whose fitness lies beyond 
Supper = (1 — r)t. We probe the tail of this fitness dis- 
tribution by again considering the deviation Gk, defined 
by Fk = 1 — Gk- This deviation obeys the convection- 
diffusion equation 



dGk 
dt 



dGk 
' dk 



D 



d 2 G k 
dk 2 



(14) 



with upward drift velocity v = (1 — r) and diffusion 
coefficient D = (1 + r)/2. The boundary condition 
G(k — vt) oc t^ 1 is set by matching the density at 
the top of the middle class with that at the bottom of 
the upper class. Consequently, the fitness distribution, 
/ = —dG/dk, follows the scaling form (Fig. 



f k (t) ~ rV 



k — vt 



(15) 



where the scaling function has the Gaussian tail 
ip(z) ~ exp(— z 2 /2), as z — > oo. 

The upper class is thus confined to a diffusive bound- 
ary layer that grows as \f~Di. From Eq. JT5J), the upper 
class contains a fraction oc X/y/t of the total population. 
Finally, for a finite population of N agents, we deduce 
from the extreme statistics criterion, Nfk ~ 1, that for 
the fittest agent fc CItrcmc ~ vt + y/tlnN. 

We comment that the deceptively simple master equa- 
tion exhibits a remarkable triple-deck structure, with a 
stationary component, followed by two transient compo- 
nents. Interestingly, the asymptotic fitness distribution 
is described by a non-analytic scaling function. 

In summary, we introduced a minimal model of social 
diversity in which the two driving mechanisms are ad- 
vancement by competition and decline by inactivity. An 
idealized but plausible social structure emerges: either a 
homogeneous society with a single lower class, or a hier- 
archical society with multiple classes. The lower class is 
always static, while the middle class and the (tiny) upper 
classes are upwardly mobile. In a hierarchical society, the 
lower and the upper classes are confined to boundary lay- 
ers that are much smaller than the dominant scale that 
characterizes the fitness of the middle class. 

There are numerous interesting questions suggested by 
this work. For example, what is the time history of an 
individual? How rigid is the social hierarchy and how 
does it depend on the population size? What happens if 



each individual is also endowed with an intrinsic fitness? 
Last, does non-trivial spatial organization emerge when 
agents move locally in space? 
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